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Abstract 

We calculate the density of states of a two dimensional electron gas located at the interface 
of a GaAlAs/GaAs heterojunction. The disorder potential which is generally created by a single 
doping layer behind a spacer, is here enhanced by the presence of a second delta doped layer 
of scatterers which can be repulsive or attractive impurities. We have calculated the density of 
states by means of the Klauder's approximation, in the presence of a magnetic field of arbitrary 
strength. At low field, either band tails or impurity bands are created by attractive potentials, 
depending on the impurity concentration. At higher field, impurity bands are created by both 
repulsive and attractive potentials. The effects of such an asymmetrical density of states on the 
transport properties in the quantum Hall effect regime, are also discussed. 
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I. INTRODUCTION 



Although the calculation of the density of states (DOS) for a two dimensional electron 



gaz (2DEG) is an old problem, it has gained interest in t 



of heterostructures grown by several groups 



re last decade because of new types 



qbsqbbbbbq. 



The heterostructures 



considered in this article are GaAlAs/GaAs heteroj unctions. They have, as usual, a delta 
doped layer behind the spacer to provide electrons in the 2D channel. Besides, they 
have another delta doped layer (Si donors or Be acceptors) inside the confining well of 
conducting electrons. Transport properties in high magnetic field have been studied on this 
kind of samples by R.J. Haug et al. jjlj. These authors observed a shift of the quantum 
Hall plateaus (for a review on the quantum Hall effect, see ref. ^j|). Merkt et al. [J| and 
Bonifacie et al. '0\ have performed cyclotron resonance experiments on these structures. 
Two main phenomena have been clearly observed. The first one concerns a new cylotron 
line, called disorder mode which is observed only for the acceptors-doped structures, 
and might be attributed to the coupling between the Landau harmonic oscillators and 
the background oscillator created by the presence of Be-atoms in the cristal. The second 
concerns a set of additional lines attributed to the energy levels created by the acceptors 
background JjJ. Kubisa and Zawadzki 3] pointed out that repulsive potentials create 
localized states when a magnetic field is present, by the combined effect of the Lorentz force 
and the electrical confinement. Although localized states created by repulsive potentials 
are at higher energy than the corresponding free Landau states, they can still trap electrons 
when the cyclotron orbit is shrinked by increasing magnetic field. These authors calculated 
the binding energies of magneto- acceptors, and showed that, contrarily to donors, acceptors 
states disappear at zero magnetic field. 

In this paper we calculate the DOS of heterojunctions having one or two delta doped 
layers. We use the Green function formalism and perform the calculations using a numerical 
method proposed earlier by Klauder[LJ|: the so-called fifth Klauder's approximation 



(KVA) |13|]. The KVA was already employed by Ando [1J] in the case of high fields, to 
obtain the DOS of a two dimensional electron gas(2DEG) in presence of a single t ype of 
impurities. More recently it was employed as well by Serre et al. [15] and Gold et a/.[l6| at 
zero magnetic field, in three and two dimensions systems. The most severe limitation of the 
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KVA is that the multiple occupancy corrections are not taken into account. Consequently 
the width of the impurity band is overestimated, resulting in a Mott density which is too 



small 
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171 ]. Another factor of imprecision is related to our choice to not incorporate 



in our model, the self energy part due to the electron-electron interaction. For sake of 
simplicity, we focused on the general asymmetrical shape of the DOS created by a given 
concentration of impurities. Several approximations for this interaction have been discussed 



in the litterature 



181 ] . but following ref. |16| . the self-energy part S ee only produces a rigid 



energy shift of the DOS. Besides, these approximations do not hold at high magnetic field, 
when the filling factor is close to one and electron-electron interaction strongly increases 
the Lande g-factor. On the other hand, we took into account the mixing between landau 
Levels (LLs) which has never been done before. Indeed, most microscopic theories Q] and 
calculations of the DOS 0, 0, 0, Q either consider only one LL, neglecting inter-LL 
mixing, or take into account only two or three LLs j0|. Consequently, the DOS is correctly 
reproduced in the strong-field limit or in the low-field limit, but is not so well known 
for intermediate or low fields, where the strength of the random fluctuations becomes 
comparable to the energy separation between LLs. 



This paper is organized as follows. In the next section, we present our theoretical treat- 
ment. We reported the results of the calculations in section III and section IV: we discuss 
the DOS for the reference sample having only one impurity layer in section III, while section 
IV is devoted to the anharmonic and asymetrical DOS. In section V, we apply this formalism 
to the transport properties in the QHE regime and we correctly reproduce the shift of the 
Hall plateaus. 



II. FORMULATION OF THE PROBLEM 



The 2D channel is considered as a homogeneous medium with a dielectric constant k and 
the charge carriers are noninteracting electrons with an effective mass m* and a charge e. 
The electrons are free in a bidimensional quantum well created by the potential U(z) of a 
heterojunction (see Fig. ^). An external magnetic field is applied in the z direction. One 
or two delta layers of impurities are introduced at a given distance from the GaAlAS / GaAs 
interface. The first doping layer has a concentration comparable to the electron concentra- 
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FIG. 1: Schematic view of the heterostructure and the fo(z) wavefunction. A delta doped layer is 
represented here, located at abscissa zq. A magnetic field is applied along the growth direction. 



tion (few 10 11 cm -2 ). The second layer of impurities, located at the distance Zq from the 
interface with a concentration iVj, create sharp potential fluctuations, either negative (for 
donors) or positive (for acceptors). Its concentration is one order of magnitude lower than 
the first one (few 10 10 cm~2). 



A. The single impurity problem 



The single impurity problem which was solved in Ref. |l2j by means of variational method, 
is here treated by a diagonalization procedure. This diagonalization method is used after- 
wards in the self consistent calculations of the Green function, for convergence purpose. 

The impurity located at point (R, 2b) creates a potential vq(t,z) = 
e 2 / (Attk,^ [z — zq) 2 + (R — r) 2 ). The three dimensional Hamiltonian of this system is 
given by 

H = -^-h + eA)\u + v , (1) 



2m* V / 

where A is the vector potential operator and p = (p±,p z ) the momentum operator. To 
solve this equation the 3 dimensional problem is first reduced to a bidimensional one 
The wavefunctions can be approximated by 

V{x,y,z)=Q(x,y)f {z) , (2) 
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where fo(z) is the wavefunction of the first subband of the quantum well with an energy Eq\ 



1 :P 2 z + U)f = E f . (3) 



2m* 

The eigenenergy E is obtained by a variational calculation and /o is approximated by the 
modified Fang- Howard trial wavefunction. Details of this method can be found in ref. Q|. 
One obtains an effective Schrodinger equation for the $ wavefunction of eigenenergy E±: 



$ = E ± $ , (4) 



where 



and 



1 r°° 

Obare(g) = pr / f (zf dz . (6) 



The screening of the bare potential by free carriers has been taken into account in ref. 



and ref. 12] by introducing the dielectric function e(q) = 1 + ^-X (q)h(q), where X is the 



Linhard function in two dimensions |2jJ and h(q) is a form factor which is due to the electron- 
electron interaction in the 2DEG. The whole problem has been solved directly by numerical 
methods. However, one should bear in mind that in the linear response approximation, the 
mathematical expression of the dielectric function when a magnetic field is applied is quite 
different from that deduced from the familiar Lindhard function in 2D. This is because both 
the energy levels and the wave functions are different. Nevertheless, we focus in this article 
on the general asymetry of the Landau levels, and not on the screening of the bare potential 
by electrons. This is why we did not use a different dielectric function. 

To calculate the DOS in the KVA approximation, we have used a diagonalization proce- 
dure instead of the variational method. Now, we apply the diagonalization method to the 
single impurity problem, in order to demonstrate first the accuracy of our procedure. We 
used the basis \N, m), N G {0, 1, ... , oo}, m 6 {— N, . . . , oo} defined by: 



(r|iV ' m) = ^ 2KP(N + \ m \)\ x e[ } Up J Up 

where I = \Jh/eB is the cyclotron length and B is the magnetic field. Because of the radial 
symmetry of the problem, the calculations can be performed for each value of m separately. 
In Fig. El we plot the calculated binding energies of magnetodonors and magnetoacceptors 



states related to the first LL. The impurities are located at the interface (zq = 0). The 
parameters are the same than those used by Kubisa and Zawadzki [12]: for a typical 
GaAs/AlGaAs heterojunction, the barrier height is 0.257 eV, the effective mass for the 
AlGaAs barrier and the GaAs well are respectively 0.073m and 0.066mo- We took a relative 
dielectric constant n r = 12.9 throughout the entire heterostructure. The electron density 
and the depletion density have been taken equal to N s = 3 x 10 11 cm~ 2 and N a = 6 x 10 11 
cm -2 respectively. The diagonalization has been performed in a basis containing a finite 
number of LLs. We have included all LL whose energy is smaller than or comparable to 
the binding energy Ej,. At B=10T, 10 LLs have been taken into account. This allowed 
us to obtain a precision better than 0.1%. At 5=0. IT, 150 LLs have been considered to 
obtain a precision of 5%. Figure 121 shows that there is no quantitative difference between 
our method and the minimization procedure used by Kubisa and Zawadzki. However the 
good convergence of the diagonalization procedure guarantees a reasonable calculation time 
for the self consistent procedure. 



The second step of our calculation procedure was to fit the true impurity potential by a 
gaussian potential which allowed us to calculate faster the binding energies, and therefore 
to optimize our algorithm in the second part of this work: 

v(r) = ^- 2 exp (-^) . (7) 

where Vq is the strength of the potential and d is its spatial extent. This correspondance 
allowed us to perform the heavy self consistent calculation of the Green function proce- 
dure, developped in the following paragraph. For this reason, we have first re-calculated 
the binding energies of Kubisa and Zawadzki by using a diagonalisation procedure with a 
fitted gaussian potential. We obtained a good fit of the binding energy of Fig. |2] by tak- 
ing Vq= 22meV and d— 10.4nm. Therefore we are in the case of short-range scatterers 



(d < 1/V2N+1, N is the LL index) [lj for a wide range of magnetic field (0< B < 6T). 
This corresponds to the position of impurities z = 0. When z = 500A, the optimized 
parameters are d = 50nm and Vq = 0.17meV. When impurities are close to the electrons, 
the gaussian potential is short range, but its becomes long range when the impurity layer is 
located at several hundreds of angstroms from the interface. 
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FIG. 2: Binding energies of magnctodonors and magnetoacceptors of the first Landau level. Pa- 
rameters are given in the text. The angular momentum m is equal to 0,1, 2 and 3. Black down 
triangles: diagonalization. White up triangles: variational calculation. At zero magnetic field, the 
binding energy vanishes for acceptor states but not for donors: the acceptor states exist only in 
presence of high enough magnetic field. 

B. Hamiltonian of the disordered system 

Now we consider the case of several impurities which are all located in the same delta-layer 
at the distance z of the hetero-interface. The Hamiltonian of the 2DEG in the presence of 
these randomly distributed impurities can be written as 



H = H ± + V, 



(8) 



with 




(9) 



and 



Vi = v(r - Ri, z Q ) . 



(10) 



Here Rj are the in-plane positions of the different impurities. 



7 



The single particle Green's function of the averaged Hamiltonian is given by the Dyson 
equation 28| : 

((G)) = 1 . (11) 

where £ is the proper self-energy part that we calculate by means of the KVA and ( ( . . .) ) 
denotes the spatial average. This approximation has been previously widely used 1J, llq 
although it is known to overestimate the width of the impurity bands [29]. This approxi- 
mation is performed into two steps. First one isolates scattering processes U occurring on a 
single impurity: 

U = Vi + ViPvi + . . . (12) 

where P = 1/(E — H±) is the Green's function of the free electron. Then, one sums over all 
the impurities: 

i=^2i i = ^2v i {i-Pv i )- 1 , (13) 

i i 

Afterwards, the bare propagator P is replaced by the dressed propagator G in the precedent 
equation: 

i = ^2v i (i-dv i )- 1 (14) 

i 

and the average value of Eq. ^2 ((• • •)) is taken over all possible positions of impurities: 

E=«t» = «X;«*(l-6i5 i )- l ». (15) 

i 

Figure represents Eq. ITol diagrammaticallv. Figure represents Eq. El Other diagrams 
are represented in Fig. |U and Fig. EJl. 

Although the problem has already been greatly simplified, the average is not known. To 
calculate the average, we follow Ando Q| who used two different basis. The first basis 
\N, X) is used for averaging over the impurities positions. The vectors of the basis are given 
by 

(r|iV, X) = ^ exp fi^ - i^j &(* - X), (16) 

where L is the length of the system, and I = \Jh/eH the radius of the cyclotron orbit. The 
functions £at(x) are written: 
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(a) x * * 




W ^ W * * /fx 



FIG. 3: Electron-impurity diagrams involved in the KVA. (a) summation of the diagrams of the 
self-energy; (b) Dyson equation; (c) examples of diagrams that are taken into consideration in the 
KVA; (d) examples of diagrams that are neglected. Dashed lines: electron-impurity interaction. 
Thin full oriented lines: bare propagators. Thick full oriented lines: dressed propagator. Crosses 
represent impurities. 



and H n (x) is the n th Hermite polynomial. 

The second basis \N, m)^ is centered on the i th impurity. This is the most convenient 
basis for summing the different scattering processes over a given impurity because of the 
conservation of the cylindrical symmetry of the impurity potential. The vectors of the basis 
are denned by 

(r|iV, m) Ri = (R ; - r\N, m) exp M ^^ ) > ( 18 ) 

where the kets \N,m) are denned by Eq. [7| The wavefunctions of the different bases are 
linked to each other by the relation 

oo 

\N,X) = J2 \N,m) n V2rt?(-l) m (R\N + m,X). (19) 

m=-N 

Finally we assume that: 

(<£|R*><H<|)>~JVi. (20) 

i 

Then, using Eq. [THl EH and EH the self-energy part can be rewritten as: 

(N,X\E\N',X') = 2Til 2 N i ^{N,m\v{\ - Gv)- x \N,m) 5 NN ,5 X x>- (21) 

m 

The two main physical parameters appear in Eq. El The first parameter is the adimensional 
concentration of impurities, c = 2irl 2 Ni which represents the number of impurities seen by 
one electron. The second parameter is the potential v. The off-diagonal elements of v are 
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responsible for the Inter-Landau-Levels Mixing (ILLM). We have already mentioned that 
the potential v can be smooth or sharp depending on the location of the doping layer. In 
the above summation, either the potential is sharp and few values of m are necessary to 
obtain a good convergence, or the potential is smooth and many values of m are needed. 
This reflects the spatial extent of the eigenfunctions which need less or more basis elements 
for their construction. 

C. The density of states 

The DOS is directly related to the averaged Green's function by the formula: 

p(E) = -13 (tr«G») , (22) 

where ((G)) is calculated by the resolution of the self-consistent set of Eq. and Eq. I2"T1 
The calculation of the DOS consists of summing matrix (whose rank is the number of LL) 
indexed by the angular momentum m over all values of the angular momentum m, in a self 
consistent procedure. In the previous section, the calculation of the binding energies was 
made by summing on the LL, for each particular value of the angular momentum m. 
The DOS is represented throughout the paper by its normalized value (NDOS) defined by 
the relation NDOS(E) = p(E)/(eB/h). NDOS(E) is such a way that: 

J NDOS(E)dE = 1 (23) 

It is represented as a function of the energy (E in mev and NDOS in mev~ l ) or as a function 
of the adimensional energy E ^ )c (the adimensional NDOS is AN DOS = NDOS x Huj c ). 

III. DENSITY OF STATES IN PRESENCE OF A SMOOTH POTENTIAL: REF- 
ERENCE SAMPLE 

We consider first the DOS of a sample having only one 5-doped layer of Si donors in the 
GaAlAs barrier, at a distance zq = —500 A from the interface. This is our reference sample: 
the disorder potential is smooth. The modification of this DOS by a strong interaction with 
a second layer of impurities (acceptors Be or donors Si) is discussed in the next section. 
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In this part, we have restricted our attention to the low coupling regime between LLs, 
when the cyclotron energy ftw c is larger than the Landau level broadening T. In this case, 
each Landau level can be treated separately. Therefore, the matrix reduces to a single real 
number, for each value of the angular momentum m. Nevertheless, a large number of m 
values is necessary to obtain the convergence of the self consistent calculation. The total 
number of m necessary to reach the convergence criterium is known, when, adding a new 
term does not change neither the shape nor the shift of the DOS. Physically, this signifies 
that the smooth potential needs a large number of states \N, m) for its correct representation. 




Energy (meV) 



FIG. 4: Density of states at B = 10 T (ground state) for a GaAlAs/GaAs heteroj unction with a 
Si — J-doped layer located at 500 A from the interface. Different curves correspond to different 
impurity concentration No = (3, 4, 5, 6).10 11 cm -2 . The surface density is Ng = 3.10 11 cm -2 . More 
than hundred "m" have been taken into account. 

Fig.4. represents the general shape of the ground Landau level whose width (overes- 
timated by our method) is increased by the presence of impurities. There are two main 
differences between this shape and the elliptic shape obtained by Ando 14, 18, 2^ or to the 
gaussian shape obtained by E. Brezin et al. j^J. First is the asymetry of the DOS: it is in 
fact enlarged on its low energy side, because of the donor caracter of impurities. Second, is 
the global shift of the DOS towards low energy. This rigid shift is due to the first term in 
the calculation of the self energy, which is generally not taken into account. The shift is a 
decreasing function of the LL index, as illustrated in Fig. |SJ This results in a anharmonicity 
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of the Landau ladder. We stress that this effect should be experimentally confirmed [3 




FIG. 5: When the first term of the self energy is taken off the equations, the Landau ladder is 
harmonic. When this term is kept, the Landau ladder is anharmonic. 

Finally, we present on Fig.|U]the DOS of the first two LL for different magnetic field values 
on an adimensional scale. We observe the disparition of the shift when the magnetic field 
is increased. The energies of the perturbed LLs are clearly approaching the unperturbed 
energies at high magnetic field: the ladder is now harmonic. 



IV. DENSITY OF STATES IN PRESENCE OF BOTH SMOOTH AND SHARP 
POTENTIALS 



1. High magnetic field 

In presence of high magnetic field, ILLM can be neglected because the adimensional 
impurity concentration is small (c = 4 x 10~ 2 for iVj=1.10 10 cm -2 at B=10T). Therefore, 
the overlap integral between two different impurities becomes negligible, because the spatial 
spreading of the localized states shrinks with increasing B. Each impurity can be considered 
as isolated from each other and the calculations shall tend to the case of one impurity 

r 

treated in ref lz 



The density of states in the high field limit is represented in Fig. [T^i and Fig. [7|3 
respectively for donors and acceptors. We observe the formation of impurity bands (IB) 
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Adimensional energy 



FIG. 6: Density of states for a GaAlAs/GaAs heteroj unctions with only one Sz-doped delta layer 
located at 500 A from the interface, for different magnetic fields. The abscissa axis is shifted by 
^r- and normalized by hu> c . The vertical axis is normalized by ^-/Huj c . 

on the low energy side in the case of donors, and on the upper energy side for acceptors. 
The energy ranges of IB are certainly overestimated due to the method used to treat the 
disorder. However, we can observe that each isolated band appearing at low density, which 
corresponds to levels of different angular momentum (m = 0, 1,2,3, ...), is exactly centered 
on the binding energy calculated by Kubisa and Zawadzki |l2|. Counting the total number 
of states of an impurity band which is entirely separated from the LL, one obtains the 
impurity density iVj. If Ebi (resp. Eb2) represents the lower (resp. higher) band edge of 
the IB, then, neglecting the spin degeneracy, one obtains: 

J p(E)dE = 2N i 

Ebi 

We observe the formation of an impurity band for densities Ni < 2.10 10 cm~ 2 . When the 
impurity concentration increases, the band width increases as well untill overlapping for a 
critical density which equals A^ c = 2.10 10 cm~ 2 at B = 10 T in our case. More precisely, 
this critical density depends on several parameters: the position of the doping layer, the 2D 
electron density Ng, and the magnetic field B. 



13 




FIG. 7: Influence of the doping density on the normalized density of states at B = 10T for the first 
LL, in the case of a GaAlAs/GaAs hetero junction having two delta doped layers. First delta-doped 
layer: concentration Np = 4.10 11 cm -2 of Si-atoms at zq = 500A Second doping layer: Si atoms 
(curve a) or Be atoms (curve b) at zq = 0. Different curves correspond to different density of the 
second layer in the range 0.1 10 10 cm- 2 < Ni < 5 10 10 cm . 

In Fig. |S1 we have reported the normalized density of states for several values of the 
magnetic field. The two sets of curves correspond to the acceptor and donor case. We notice 
that the global shift represented in Fig. |H] is not strongly modified by the additional low- 
density impurity layers. Then, increasing the magnetic field strength, leads to the formation 
of a separated impurity band for the N = LL. At low field, this impurity band collapses to 
form a band tail. For the first excited LL, no separated impurity band are created, because 
of too small binding energies. 

2. Intermediate magnetic field 

We comment in this paragraph, the overlap of LLs as a consequence of their cou- 
pling when the parameter c overpasses 1 (by increasing N). The DOS calculated for 
B = 0.5 T, is represented in the case of donors in Fig. Eland in the case of acceptors in Fig. [Till 

In both cases the DOS is enlarged in the low energy side by the long range potential 
fluctuations of attractive Si parents donors at 500 A from the interface. The addition of 
short range potential fluctuations at the interface enlarges the DOS differently in case of 
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FIG. 8: Density of states for different magnetic field values, in the case of a GaAlAs/GaAs hetero- 
junction having two delta doped layers. First delta layer: Si atoms at zq = 500A. Second doping 
layer: a) Si atoms at z$ = ; b) Be atoms at z$ = 0. 

acceptors or donors. Attrative Si impurities enlarge the DOS on the same side than the 
parent donors, whereas Be impurities enlarge the DOS on the other side. 

Up to the critical value c = 1, the degeneracy of the first LL is large enough to contain 
all impurities states. When the density of the impurity layer increases ( c > 1), the LLs 
overlap as shown in Fig. The first LL remains separated from the higher LLs for a wider 
range of concentration values, because the shift is more important for low LL index N, as 
previously remarked. 

Fig. is dedicated to the acceptor case. Now, the first LL remains separated from 
the Landau ladder even at the higher density of Aj = 8.10~ 10 cm -2 and the merging of LLs 
begins by the high values of N. This difference between acceptors and donors is related to 
the difference which is observed at B = 0T : there only exists donors states at very low 
B. Moreover, in the case of attractive potentials, the coupling between LLs increases the 
binding energies of the lowest LLs and decreases the binding energy of the highest LLs, 
resulting in a donor state of finite binding energy at B — 0T (see Fig. 2). 
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FIG. 9: Donors case: density of states for a GaAlAs/GaAs heteroj unctions having a delta layer of 
Si at 500 A from the interface, and another delta layer (density Ni) of Si located at the interface. 
The magnetic field is B = 0.5 T. 



3. Low magnetic field 

Fig. ITT1 shows the DOS in the case of donors impurities at B = 0.1T as a function of 
energy. Different values of the donor concentration Ni have been considered, and the rigid 
energy shift of the different curves has been removed. The energy reference is fixed at the 
center of the first unperturbed Landau level. The high energy side of some curves has been 
cut off for clarity. In the acceptor case, which is not represented here, no IB below the 
conduction band appear . 

At low field, ILLM becomes dominant: LLs overlap and tend to the flat Conduction 
Band (CB) of a 2DEG at B — 0T. At such a low B, the degeneracy of one LL is 
smaller than 2Ni and several LLs contribute to the formation of the IB (ILLM). Fig. ITTIrep- 
resents the donor case. In agreement with Fig. 2, only one IB exists, corresponding to m — 0. 
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FIG. 10: Acceptors case: density of states for a GaAlAs/GaAs heteroj unctions having a delta 
layer of Si at 500 A from the interface, and another delta layer (density iVj) of Be located at the 
interface. The magnetic field is B = 0.5 T. 

The critical density (c = 1) is represented by trace 7. Above this value, the IB and the 
CB merge and form a band tail. At vanishing concentration c, the position of the IB is again 
in perfect agreement with the binding energies calculated in section 2 (see Fig. |2J). 

These results are in good agreements with the theoretical results of Gold, Serre and 
Ghazali 16] who calculated the DOS of a 2DEG in the case B = 0T within the same 
Klauder's approximation. 

We do not observe the formation of a second IB that could originate from states of 
higher momenta. As a matter of fact, our potential does not allow the formation of excited 
impurity bands at vanishing magnetic field because all the binding energies vanish at 
B = 0T, except for m = 0. The disappearance of the IB occurs at a concentration N£ = 0.5 
10 10 cm~ 2 . Gold et al. pointed out that the enhanced Thomas-Fermi screening they used 
implied large impurity bandwith. However the critical metal-insulator densities we found 
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FIG. 11: The density of states in the case of magnetodonors, at B = 0.1T, for different impurity 
concentrations. The curves labelled from 1 to 11 correspond to concentrations 0.01, 0.025, 0.05, 
0.1, 0.15, 0.2, 0.25, 0.5, 1, 2.5 and 5.0 x 10 10 cm~ 2 respectively. Curve 7 corresponds to the 
concentration Ni = 0.25 10 10 cm -2 , for which c = 1. Inset: enlargement of the IB at low donor 
concentrations. The flat conduction band degeneracy is equal to -^p. 

are so small that they are more probably related to the errors induced by the Klauder V 
approximation itself |l7|. 



V. QUANTUM HALL EFFECT 

It has been established jlj in GaAs/GaAlAs heterostructures that 2DEG perturbed by 
acceptors or donor impurities exhibits shifts of the quantum Hall plateaus relative to the 
line for the classical Hall resistance. If acceptors are added, then the plateaus shift to lower 
filling factor v. At the contrary, if donors are added to the 2DEG, then the plateaus shift 
to higher v. These shifts could be succesfully explained by the asymmetry of the DOS in 
each LL. Because our theoretical approach takes into account ILLM, the Hall resistivity Rh 
at both high and low magnetic fields can be calculated and therefore a direct comparison 
between the classical Hall effect and the shifts of the plateaus is possible. 
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FIG. 12: Three curves corresponding to three different "samples". Reference is represented by a 
full line. The Hall resistance of the acceptor-doped sample is shifted towards low filling factors 
(higher B), whereas it is shifted towards high filling factors (lower B) in the case of donors. 

For sake of simplicity we assume that only one state per LL is delocalized whose energy E n 
is given by the maximum of the n-th LL j^l| . Furthermore we assume that this state exists 
even in the presence of strong disorder and ILLM and that it contributes to the conductivity 
with e/h. Within these approximations, the adiabatic conductivity is given by 

n 

where / is the Fermi distribution and Ep the Fermi energy. Figure fHZl shows the calculated 
Hall resistance as a function of the magnetic field B. For all of the three cases of Figure fT2"l 
a 5-doped layer located at 500A from the interface, with a concentration Np = 4.10 11 cm -2 , 
contains the parent donors. A second 5 layer of donors (acceptors) is added at zq = 0, 
and shifts the curve towards low (high) magnetic fields, compared to the reference case. 
The figure reproduces well the shift of the Hall plateaus, and must be compared to Fig. 1 
and Fig. 2 of Ref. [jj. Experimentally, the shift of the plateaus is systematically much less 
pronounced for donors than for acceptors, and this effect is reproduced here. The reason 
lies in the fact that the remote layer is always doped with donors. At low field the three 
curves of Fig. H~2l merge into the same classical line. 
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VI. CONCLUSION 



We have used the multiple scattering approach proposed by Klauder and the averaging 
procedure proposed by Ando to calculate the density of states of a disordered two dimen- 
sional electron gas in the presence of a magnetic field of arbitrary strength. The impurities 
are randomly distributed in a layer at a given distance from the electron gas. Each impurity 
creates a Coulomb potential. The screening of these potentials due to the 2D electrons was 
taken into account. The impurities were either donors or acceptors. 

At low field, in the donor case, the density of states exhibits either an impurity band or 
a band tail, depending on the concentration. In the acceptor case, there is neither impurity 
band nor band tail for low field. 

At higher field, one observes that impurity bands split from the Landau Levels. The 
situation between donors and acceptors becomes symmetric only at high fields, where the 
mixing between the Landau levels can be neglected. 

Our results are in qualitative agreement with transport experiments on magnetic-field- 
induced metal-nonmetal transition in GaAs-GaAlAs heterostructures 3^|, with recent 
magneto-photoluminescence experiments |33] and with transport experiments in the Quan- 
tum Hall regime [lj]. This model can be extended to calculate magneto-transport and 
magneto-optical properties of the 2D electron gas. 
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